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We study the enhancement of the ferromagnetic relaxation rate in thin films due to the adjacent 
normal metal layers. Using linear response theory, we derive the dissipative torque produced by the 
s — d exchange interaction at the ferromagnet-normal metal interface. For a slow precession, the 
enhancement of Gilbert damping constant is proportional to the square of the s — d exchange constant 
times the zero-frequency limit of the frequency derivative of the local dynamic spin susceptibility 
of the normal metal at the interface. Electron-electron interactions increase the relaxation rate by 
the Stoner factor squared. We attribute the large anisotropic enhancements of the relaxation rate 
observed recently in multilayers containing palladium to this mechanism. For free electrons, the 
present theory compares favorably with recent spin-pumping result of Tserkovnyak et al. [Phys. 
Rev. Lett. 88,117601 (2002)]. 



I. INTRODUCTION 

Ferromagnetic multilayers have attracted much atten- 
tion recently because of their applications in spintron- 
ics and high density magnetic recording devices. The 
present paper is concerned with magnetic relaxation in a 
ferromagnetic film (F) imbedded between nonmagnetic 
metallic layers {N). In particular, we study the enhance- 
ment of the Gilbert damping in N/F/N sandwiches as 
compared with that of a single ferromagnetic film. The 
Gilbert damping constant G is defined by the Landau- 
Lifshitz- Gilbert (LLG) equation of motion [1,2] 
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where M is the magnetization vector, Ms is the satura- 
tion magnetization, and H^f / is the effective field which 
is given by Heff = —dE/dM where E is the Gibbs free 
energy. The gyromagnetic ratio 7 is a negative quantity 
— g/i where g is the spectroscopic splitting factor, and 
/is is the Bohr magneton. The second term on the R. H. 
S. of Eq.(l) represents the dissipative torque, the mag- 
nitude of which is proportional to G. The LLG equation 
describes well both the static and the dynamic properties 
of ultrathin ferromagnetic films and the N/F/N sand- 
wiches (for a review see Ref.[3]). The experimental and 
theoretical aspects of spin relaxation in multilayers are 
covered in Ref.[4]. 

Recent studies of the FMR linewidth for ultrathin films 
[5,6] show that the constant G is enhanced when a non- 
magnetic metal is deposited on the ferromagnetic film. 
This effect has been predicted in 1996 by Berger [7] in his 
study of spin transfer between the polarized conduction 
electrons and the ferromagnetic film. Berger shows that 
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there is an enhanced electron-magnon scattering caused 
by isotropic s — d exchange taking place at the F/N in- 
terface. Being a surface effect, the enhancement of the 
Gilbert constant is inversely proportional to the thick- 
ness of the ferromagnetic film. This unique feature is 
observed in recent measurements of the additional FMR 
linewith on permalloy-normal metal sandwiches [8] , and 
on the double ferromagnetic layer structures as well [9]. 
It should be mentioned that already in 1987 a study of 
FMR linewidth showed an appreciable increase in the 
Gilbert damping with decreasing thickness of the Fe ul- 
trathin films grown on bulk Ag(OOl) substrates [10]. 

A novel mechanism for additional Gilbert damping 
in N/F/N structures has been recently proposed by 
Tserkovnyak et al.[ll]. These authors calculate the spin 
current pumped through the N — F contact by the pre- 
cession of the magnetization vector M{t). The theory is 
based on extending the scattering approach of parametric 
charge pumping by Brouwer [12] to spin pumping. 

Like the theory of Berger [7], the additional damping 
of Rcf.[ll] scales inversely with the thickness of the fer- 
romagnetic film indicating that only the F/N interface 
is involved. However, the expression for excess G, which 
we call G", differs considerably from that of Ref.[7]. In 
particular, G' vanishes with vanishing s — d exchange 
splitting. An attractive feature of this theory is that it 
links G' to the transport properties of the interface. Due 
to exchange polarization of the F/N contact, the reflec- 
tion (r) and transmission (t) coefficients at the interface 
depend on the orientation of the conduction electron spin 
with respect to the magnetization direction of the ferro- 
magnet. The formula for G' involves differences such as 
Ar = r'f —rK Interestingly, similar quantities play a role 
in the theory of interlayer magnetic coupling by Bruno 
[13]. For instance, the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) coupling [14] between two-dimensional layers 
can be obtained by calculating the interlayer coupling 
energy in terms of the reflection coefficients and 
of the layers. These coefficients are obtained by solving 
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a simple problem of scattering by J-function potential. 
In the limit of weak exchange splitting (compared to the 
Fermi energy) , the derived intcrlayer coupling agrees with 
the RKKY result of Yafet [15]. 

These considerations prompt us to take another look at 
the theory of enhanced relaxation in multilayer. The fact 
that the RKKY theory [13] involves transport properties 
at the interface similar to the spin pump theory of G' [11] 
suggests that a suitable generalization of RKKY theory 
to time-dependent magnetization M{t) may unravel the 
needed dissipative torque of Eq.(l). We note that the 
standard approach to RKKY coupling is to use a lin- 
ear response theory [14,16], and calculate the conduction 
electron spin density induced by the contact exchange 
potential. In applications to interlayer coupling between 
ferromagnetic layers with time-independent magnetiza- 
tion vectors, it is the static spin susceptibility of the elec- 
tron gas which determined the coupling. In the present 
paper, we consider a response to a slowly varying time- 
dependent s — d exchange potential. 

Owing to the dissipative part of the spin susceptibility, 
the spin density induced by the precession of M{t) will 
have a component that is out of phase with M{t). Hence, 
such a component will have a time-dependence given by 
dM/dt. The induced spin density has both local and 
nonlocal consequences for the dynamics of multilayers. 

First, let us focus on the local effect in the N/F/N sys- 
tem. It is instructive to invoke the analogy to radiation 
damping of charged particle in classical electrodynamics 
[17]. Thus, we put the dissipative torque of Eq.(l) in 
correspondence with the radiation reaction force acting 
on the particle. The first term of Eq.(l) corresponds to 
the external force. In view of this analogy, we rewrite the 
dissipative torque in terms of a reaction field Hj. 
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where H^!'"' represents the dissipative part of Hj.. For 
the N/F/N system shown in Fig.l., we find a reaction 
field (see Appendix A) 



(3) 



where J is the s — d exchange coupling constant, a 
is of order of the lattice constant, and d is the width 
of the ferromagnetic film. The quantity {s{x = 0,t)) 
is the spin density induced at the F/N interface by the 
s — d exchange interaction (x being the distance from 
the interface). The expression (3) is quite general in the 
sense that it can be used for both ballistic and diffusive 
cases. In the present paper we confine ourselves to the 
ballistic case and calculate the induced spin density using 
linear response theory [16]. 

Assuming slow precession of M{t), we find two contri- 
butions: One that is proportional to M{t) with a coeffi- 
cient given by the real part of the local spin susceptibility 



at zero frequency. If the spin susceptibility is anisotropic, 
this term leads to an anisotropic shift of the FMR fre- 
quency. For isotropic susceptibility, the quantity on the 
R.H.S. of Eq.(3) is a vector parallel to M{t) and the 
corresponding torque (2) vanishes. 

The second contribution to {s(x = 0, t)) is proportional 
to dM/dt. The coefficient of proportionality is the fre- 
quency derivative of the imaginary part of the suscepti- 
bility taken at = 0. Like the real part, this quantity is 
generally anisotropic. According to Eq.(2), it produces a 
dissipative torque leading to an anisotropic G'. Explicit 
evaluation of this term for an isotropic nonintcracting 
electron gas yields a formula for G" which compares fa- 
vorably with the spin pump theory of Tserkovnyak et 
ah [11]. 

The proposed formulation allows us to incorporate in- 
teractions between the electrons in the N region. The 
generalized Hartree-Fock approximation shows that the 
q, uj dependent spin susceptibility of the interacting elec- 
tron gas is enhanced compared to the free electron case 
[16]. Using these results, we find that the part of 
that contributes to the FMR frequency shift is enhanced 
by the Stoner factor 



Se = [1- UNiep)]- 



(4) 



where U is the screened intraatomic Coulomb interac- 
tion and N{eF) is the electron density of states, per atom, 
at the Fermi level [16]. On the other hand, we find that 
G' is enhanced by a factor of ^l;. The Stoner enhance- 
ment is thought to be large in metals such as palladium 
and its alloys. Recent results using multilayers with Pd 
layers as a spacer show a significant enhancement of inter- 
face damping exhibiting a fourfold anisotropy in keeping 
with the present theory [18]. 

Let us turn to the nonlocal consequences of the dynam- 
ically induced spin density. This is relevant in the case 
of two time-dependent magnetizations Mi{t) and M2{t) 
separated by a nonmagnetic spacer. Similar to the local 
effect, the magnetization Mi{t) at xi induces at the po- 
sition X2 of M2{t) a spin density oscillation with both 
the in phase (oc Mi{t)) and the out of phase component 
(cx: dMi/dt). The interaction of the in-phase component 
with M2{t) yields the standard interlayer RKKY cou- 
pling [15]. The out of phase component leads to a dissipa- 
tive interaction of the form dM\{t)/dt.M2{t)f{x2 — x\). 
A similar form, with the indices 1 and 2 interchanged, 
follows by taking into account the spin density induced 
at X\ by precession of M.2{t)- Such terms may be of im- 
portance in FMR studies involving dynamically coupled 
magnetic films. However, these topics are beyond the 
scope of the present investigation which focuses on dy- 
namics of single ferromagnetic film in contact with non- 
magnetic metal. 

The paper is organized as follows. In Sec. II we de- 
rive the spin density induced by a two-dimensional layer 
of precessing spins, and the corresponding reaction field. 
Sec. HI focuses on the dissipative part of the reaction field, 
and the excess damping G' for an isotropic electron gas. 
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An expression for G' in a free electron model is derived 
for N/F/N structure with N layers of both infinite and 
finite thickness. The enhancement of G" due to electron- 
electron interactions is considered within the generalized 
Hartree-Fock approximation. In Sec. IV we establish a re- 
lation between the spin pumping theory of Ref.[ll] and 
our free electron result for G". The role of spin sink in 
theories of Gilbert damping is discussed in Sec.V. The 
reaction field for the N/F/N structure is derived in Ap- 
pendix A. The role of spin relaxation in the theory of G' 
is considered in Appendix B. 



II. DYNAMIC RKKY 

We consider a two-dimensional layer of aligned spins 
imbedded in normal metal. Notice that a similar model 
was used by Yafet [15] to calculate the interlayer coupling 
for time-independent magnetizations. Here we take into 
account the time-dependence of the precessing magneti- 
zation. 

Our task is to derive the conduction electron spin den- 
sity induced by the s—d exchange interaction taking place 
at the magnetic layer. The Hamiltonian of the conduc- 
tion electrons is 



sheet 

H = Ho-JY^ / d^rS'^'\t).s{r)6'^{r - n) (5) 

i 

where Hq represents the Hamiltonian of the electrons in 
the absence of the ferromagnetic layer, and the second 
term is the s — d exchange interaction with the layer. 

J is the s — d exchange coupling constant, 5''*'(t) is 
the classical spin at the location r-i, and s(r) is the spin 
density operator for the electrons in normal metal. 

For a uniform precession of aligned spins, we have 



S^^{t) = -M{t) 
7 



(6) 



where fl is the volume of the unit cell. 

The second term of Eq. (5) acts as a time-dependent 
perturbation that will be treated using the linear re- 
sponse theory [16]. Thus, the expectation value of the 
/i-component of the induced spin density is 

sheet n p OO 

{s^{r,t)) = JY. / dt'x^^{rt,r't') 

xSj^\t')S^{r'-ri) (7) 

where 

X^Art,r't') = ^e{t-t'){[s^{r,t),s4r',t')]) (8) 

is the retarded spin correlation function (susceptibil- 
ity)[16], Q{t — t') being the unit step function. 



Now, we evaluate the R.H.S. of Eq.(7) assuming a slow 
precession. This approximation is valid as long as the 
precessional frequency is small compared to the relevant 
excitation frequencies of the conduction electrons. We 
note that the spin pump theory requires a similar condi- 
tion to hold [11, 12]. 

Performing the r'-integration, and applying the com- 
mutation law for the convolution, Eq.(7) reads 

sheet „oo 

(s^(r,t)) = J^ / dt'S\:\t-t')x^,{r,ri,t') (9) 
i 

For slow precession, we write 



S^\t-t')^S\:\t)-t'^-^ (10) 

Introducing this expansion into Eq.(9), we obtain with 
use of Eq.(6) 



(s^(r,i)> ~ — hm 



sheet 



dt 
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where 



/oo 
rfte'"*X/..(r-,ri,t) (12) 
-oo 

We note that the second term on the R.H.S. of Eq.(ll) 
follows by letting the derivative of the real part of the 
susceptibility with respect to to equal to zero at w = 0. 
This is consistent with the reality condition implying that 
that the real part of the susceptibility is an even function 
of io. The first term on the R.H.S. of Eq.(ll) corresponds 
to the static RKKY result of Yafet [15]. The second 
term is the dynamic generalization of RKKY for a slowly 
precessing magnetization. 

For an infinite medium, the sheet sums in eq. (11) 
can be evaluated by invoking translational invariance and 
Fourier transforming XiJ.v{'r,ri,Lo) to q-space [13]. Thus, 
we define a generic sum 



sheet 



d'^J 



X^^,{r,Lo) = Xi.Ar,ri,Lo) - j ^^^^^ 

sheet 

^ exp[iq.{r -ri)]x^^{q,Lo) (13) 



from which both terms of Eq.(ll) can be deduced. 
To perform the g-integral on the R.H.S. of this equa- 
tion, we set 
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9 = 9ii +g± 
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where q|| is confined to the sheet, and gj^ is perpen- 
dicular to the sheet. For a square sheet of area L^, we 
have 

f d^q _ f d'^QW f dq± _ 1 /" ^g-L f.^s 
J (27r)3 J {2nY J 27r 2^ J 2tz ^ > 



Assuming a continuous distribution of spins, the sheet 
sum in Eq. (13) is given by 

sheet 

^ exp[ - i(q|| + qj_).ri] = N^dq^^fi (16) 



and the dissipative part of the reaction field Hf^ is ac- 
cording to Eq.(19) given by 



2J^Qans 



X lim 



9 r dq^ . ; 



72 fi 
dM{t) 



(20) 



dt 



Introducing this result in the R.H.S. of Eq.(2), we ob- 
tain the damping enhancement constant G' 



G'~2j20n,M2(-) lim 

a ui—>o 



d_ 
dto 



dq 
2^ 



Imx{q,oj) 



(21) 



A. Independent electrons 



where Ns is number of spins in the sheet. Using Eqs. 
(14-16), the R.H.S. of Eq.(13) is evaluated with the result 



= n. 



dq± 
27r 



exp(iqr_L.r)x^,.(q_L,w) (17) 



where Ug = Ng/L'^ is the sheet density. Owing to the 
symmetry of the model, the quantity X^^{r,uj) depends 
only on the distance x from the sheet. 

Also, the induced spin density is a function of x. Using 
Eqs.(ll), (13) and (17), we obtain 



(s^(a;,i)) = — lim X^^{x,ijj)M^{t) 
dImX^^{x,uj) dM^{t) 



(18) 



du) 



dt 



First, we evaluate the expression (21) by disregarding 
the electron- electron interaction in the N regions. How- 
ever, a finite splitting A of the t and [ spin bands is 
assumed. The external magnetic field of the FMR ex- 
periment is one source of this splitting. For a system of 
infinite size, this splitting establishes a lower cutoff on 
the wave vector q. As shown below, this cutoff is es- 
sential to prevent the logarithmic divergence of Eq.(21). 
Due to the spin splitting, the susceptibility develops some 
anisotropy. Since transverse components of the reaction 
field (19) contribute to the dissipative torque, we need to 
use in Eq.(21) the transverse susceptibility [16] 



(0)/ \ 
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+ A 



(22) 
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Using this result and Eq.(17) in Eq.(3), the /x- 
component of the reaction field is 



2j2flans 
Hr,i^{t) ~ — — hm 

j^d w^o 

d 



dq 
2^ 



X^.(g,co)M,(t) (19) 



dq^ , .dM^{t) 



The first term on the R.H.S of this equation contributes 
to the torque [M x Heff] only if X/^i^ is anisotropic. In 
this case an anisotropic shift of the FMR frequency en- 
sues. On the other hand, the second term contributes 
to the dissipative torque that is non vanishing for both 
isotropic and anisotropic susceptibility. 



III. GILBERT DAMPING 

In what follows, we consider the dissipative torque for 
an isotropic electron gas. Thus, Xfi.i^iQj'^) = x(9) 



Expanding the Fermi functions, we have 



(24) 



Using Eqs. (23) and (24), the imaginary part of Eq.(22) 
becomes 



t2 /'OO 



Imx^\q,oj)':^ — I dkk^-^[ dOsine (25) 
167r Jo dek Jo 

X Ah |(5( fej- A 



m 



m 



where we used polar coordinates to perform the k- 
integration. Performing first the integration over polar 
angle 6, we have 



d^sin^l A+^^V^/u^-A-^') (26) 



m 



hkq 



e fe- 



rn 
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Introducing this result into Eq. (25) , and converting the 
fc-integral to e^-integration, we obtain at T = 



xe 



2me 



e k 



\hLu — A\m 
\hLu — A\m 



The unit step function on the R.H.S. of this equation 
equals one for q > qi, and zero for q < q\ where 



qi 



\hijj — A|m 



(28) 



Thus, qi acts as a lower cutoff in the ^-integral of 
Imx{q,oj). Since \hu; — A| -C ep, the upper cutoff is 
given by q2 « ikp- We then get using Eq.(27) 



dq ^ (0) . X 



In 
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(29) 



IGn^h \tko - A| 



Introducing this result into Eq.(21), we obtain 



G' 



(30) 



where we assumed a cubic lattice to write Sluga = a^. 
Strictly speaking, G" in this equation is not a simple 
Gilbert damping since its strength depends on the ap- 
plied magnetic field, but it is a weak dependence. 

When A vanishes, G' exhibits a logarithmic diver- 
gence. This is due to the fact that precessing magne- 
tization radiates into the A'^-region a very large number 
of electron-hole pairs of low energy. This is reminiscent 
of the infrared divergence in quantum electrodynamics 
where a large number of low energy photons is radiated. 
Endowing the photon with a small mass cures this diver- 
gence [19]. In an analogous way, a nonzero spin splitting 
of the electron energy levels suppresses the mimber of 
low energy electron-hole pairs that are generated by the 
precessing magnetization. 

For A ~ gi-b/sH, we have 



In 
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We see that it is the magnetic length 



obtain from Eq.(32) Ih — 4x lO^a. In real samples, the 
layers are usually considerably thinner than this length. 
As shown in the next subsection, the relevant cutoff is 
then determined by the boundary conditions (b.c.) at 
the outer surfaces of the sandwich. 

It is interesting that spin relaxation of the conduction 
electrons can provide an effective gap and a finite cut- 
off length even in the limit of infinitely thick N layers. 
As shown in Appendix B, the resulting cutoff length is 
qi^ K vpTs where Ts is the spin-relaxation time. The 
corresponding value of the effective gap is Ae// « h/Tg. 
The same quantity determines the linewidth of the elec- 
tron paramagnetic resonance of conduction electrons [23] . 
Its magnitude is presumably not negligible in comparison 
with the magnetic splitting qhbH. 



B. Finite normal layers 

We consider a ferromagnetic sheet imbedded between 

two A'' layers each of thickness D. Strictly speaking, in a 
finite system, the translational invariance invoked in the 
evaluation of the sheet sum (13) is broken. This compli- 
cates the calculation of the induced spin density. Never- 
theless, a simplified approximate evaluation of Xfj,,y{x,u)) 
can be carried out if x ~ 0. In this case, the transla- 
tional invariance is restored locally since the boundary 
at X = D plays a small role. The b.c. at a; = is taken 
in the form 



(33) 



This condition follows from assuming that there is an in- 
finite potential step at the normal metal- vacuum bound- 
ary. Thus, the A^-electron wave functions are forced to 
zero at the boundary and so is the magnetization den- 
sity. Applying this b.c. to Eq.(18), the finite- Z) version of 
Eq.(17) reads 



X^^{X,U!) - ^ X! ^^^(inX)Xt^iyiqn, 



(34) 



n=l 



where (/„ = 7r(n — l/2)/£>, and rimax — D/a since it is 
the lattice spacing which determines the highest value of 
qn- In view of this, Eq.(21) is changed to 



G' ~ ^ j2^„,m2 lim 
dD " 



d_ 
dco 



(35) 



For noninteracting electrons without spin splitting, we 
have 



he 



(32) 



which plays the role of q^^ ^ in case of infinitely wide A'^ 
layers. Taking H = 10^ gauss and putting kp — ^/a 



we 



16ttH 



(36) 



For D/a » 1, we use the definition of the Euler constant 
7b ~ 1.78 to obtain 
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TT 



(37) 



Using Eqs. (35-37), we obtain 



(38) 



As expected, the boundary conditions in a finite slab 
imply a cutoff qi pa D~^. 

Wc now make an order of magnitude estimate of 
Eq.(38) for an iron film of thickness d = D — 10a where 
a = 4 X 10~®cm. The constant J can be estimated by 
relating it to the atomic exchange integral Jgd 



J: 



O.leV, and Ms 
-1 



(39) 

IJxlO^gauss, Eq.(38) 



Taking J^, 

yields G' « 10**s~^. This agrees with the interface damp- 
ing observed recently in the double layer structure by 
Urban et al.[9]. It should be pointed out that the sec- 
ond ferromagnetic layer in this experiment plays a crucial 
role in establishing the spin sink needed to prevent spin 
accumulation in the N layers( see Sec.V.) 



C. Electron-electron interactions 

We now calculate G' by taking into account interac- 
tions between electrons in the normal metal. The gen- 
eralized Hartree-Fock approximation for the Hubbard 
model yields the following expression for the transverse 
susceptibility [20] 



(40) 



where U = A^U/h?, and U is the screened intraatomic 
Coulomb energy. Using this formula, we have 



lim 



d 



*o duj 



ImxT{q,uj) 



(9,0) 



(41) 



X lim 



d 



ImXT\q,<^) 



To simplify the evaluation of the g-integral of this 
quantity, we take advantage of the weak dependence of 
the static susceptibility on g for g < qp, and make the 
approximation 



X^°^(9,0)«xV"^(0,0) = — iV(e^) 



(42) 



where N{ep) is the density of states, per atom, at the 
Fermi energy. Using Eqs. (41) and (42) in Eq.(35), we 



obtain the enhancement, G', for interacting electrons in 
finite N layers 



G" 



{JMsamSE? 
87r2M 



ln(i:>/a) 



(43) 



where Se is the Stoner factor defined in Eq.(4). For pal- 
ladium, we have Se ~ 10. Thus, large values of G' are 
expected for sandwiches containing Pd as normal layers. 
Mizukami et al.[8] measured the Gilbert damping con- 
stant in N/F/N sandwiches with F being a thing film of 
permalloy (Py). These measurements show that G" for 
Pd/ Py/ Pd system is well above that for the Cu/ Py/Cu 
system. However, it is about twice as smaller than G' 
for Pt/Py/Pt. This seems to contradict our Eq.(43). 
We attribute this disagreement to spin accumulation in 
the normal layers. Since the spin-orbit coupling constant 
of Pt is larger (by a factor of 3) than that for Pd, the 
spin lattice relaxation rate in Pt is an order of magnitude 
stronger than that in Pd. As pointed out by Tserkovnyak 
et al [11], the spin accmnulation takes place when the 
spin relaxation rate is small. The fact that experimen- 
tal ration Gp^/G'p^ is less than 10 indicates that Stoner 
enhancement in Pd is not excluded. 

More convincing evidence of Stoner enhancement in 
interface damping comes from recent FMR studies on 
20Au/40Fe/40Au/5Pd/[Fe/Pd]5/UFe/GaAs{001) 
samples [18]. Compared to a single layer structure, these 
samples show a G' that is enhanced by a factor of four 
and exhibits a strong fourfold in-plane anisotropy. Ap- 
parently, the presence of a second Fe layer provides an 
efficient spin sink. Thus G' is determined by the exchange 
enhanced susceptibility rather than the bottleneck due to 
a weak spin-lattice relaxation in the Af-layers. 

We now digress, for a moment, to consider the en- 
hancement of the FMR frequency shift due to interac- 
tions. Applying Eq.(40) to the static anisotropic sus- 
ceptibility, and making the approximation (42) , we have 
Xfii^il,^) ~ 'S'£;x^°2(g, 0). Using this result in Eq.(19), we 
see that the frequency shift for the interacting electrons 
is Se times that for the independent electrons. This pre- 
diction could be used, in conjunction with the data for 
the anisotropic G' to clarify experimentally the role of 
interactions in the FMR of multilayers. 



IV. RELATION TO SPIN-PUMPING THEORY 

We now show that for free electrons there is a similarity 
between our formula (30) for G" and the spin-pumping 
theory of Tserkovnyak et al. [11]. According to these 
authors, the excess damping produced by pumping of 
spins into adjacent N-layers is G' = 'yM^a' where 



a 



AnMsL'^d 



(44) 



where gj^ is the Landc factor, fig is the Bohr magneton. 



and A. 



(L) a(R) 



are the interface parameters for the left, 
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right A'^-Layers, respectively. In terms of the elements 
of the 2x2 scattering matrix, for a symmetric N/F/N 
sandwich, these parameters are 



4 



2 / J II mn 

mn 



'T 



(45) 



over the in-planc wave vectors fe|[. Converting the fey 
sum to a two-dimensional integral, we have 



A 



dk\\k\\ 



(52) 



Evaluating this integral, we get Ar ~ {L^P"^ /'!t)F{I3) 
where 



where (rj,^„, r^^) and (tj^mi ^mn) ^re the reflection and 
transmission coefficients for electrons with up and down 
spins. The expression (45) is to be evaluated with the 
transverse modes (m, n) taken at the Fermi energy. 

Following Bruno [13], we consider the scattering of A''- 
elcctrons by a ferromagnetic monolayer. Due to the con- 
servation of transverse momentum, rmn = fmSmn where 
fm = ro{k±), k± being the component of the electron 
wave vector perpendicular to the monolayer. 

The reflection coefficients ro{k±) are found by solving 
the one-dimensional scattering problem for the potential 



v{x) = VoS{x) 



(46) 



where Vq is given by the interface coupling constant J, 
and by the magnitude of the atomic spin S of the ferro- 
magnet 



vo = ±^JSns 



(47) 



where the (-I-, — ) signs correspond to the (J,, j) electron 
spins, respectively. 

The reflection coefficients for this problem are [13] 



F(/3)~ln.^^+^' 



/32 kl + /32 



(53) 



Inserting Eqs.(52) and (53) into Eq.(44) and expressing 
with use of Eqs.(47) and (48), we obtain 



7(mgJn,)2 



(54) 



We can bring this result closer to the form of Eq.(30) by 
letting Hs = a~^, and Mg = ^Sa~^ yielding 



2^2 



Mia 
7 



(55) 



Furthermore, assuming /3 <C kp, the function F{P) can 
be approximated by 



i^(/3)«21n-^«21n4^1 
1.65/3 Jsd 



(56) 



^From Eqs. (54-56) we have 



-i/3 



k±+iP 
k±-i(3 



where 



/3 = 



mvo 



The transmission coefficients are [20] 

t ik± 



^ ik±-l3 
A _ 

'~ik^+l3 



Eqs. (48) and (50) imply 



4/32 (fc 



(k 



(48) 



(49) 



(50) 



(51) 



where we applied the identity (fcjj + k'j_)^p 
this result, the transverse mode sum (45) becomes a sum 



kp. Using 



G' = -iM,a' 



This equation shows a remarkable similarity with the ex- 
pression (30). Note, however, that the logarithmic terms 
do not match. If we consider an infinite system, and ig- 
nore the cutoff due to the magnetic length (32), the gap 
A vanishes , and Eq.(30) becomes logarithmically diver- 
gent. On the other hand, Eq.(57) shows that there is 
an effective gap, A « Jgd corresponding to a finite cut- 
off qi ~ kpiJsd/ ^f)- Thus, the spin-pumping theory is 
infrared-divergence free. For real finite size systems, this 
difference is only of academic interest since the cutoff 
produced by the boundary conditions, qi « -D^^, yields 
a logarithmic term that is of the same order of magnitude 
as that in eq.(57). 

The presence of an effective gap, A « Jsd, in the spin 
pumping theory is presumably linked to the fact that, in 
contrast to linear response theory, it is of infinite order 
in the coupling constant J. This is seen in Eq.(54) where 
F{P) is a nonlinear function of /? given in Eq.(53). This 
kind of nonlinearity is generic in the scattering approach 
to transport(see Eq.(51)). In fact, Bruno [13] derives an 
exact expression for the static RKKY coupling that goes 
beyond the hnear response result of Yafet [15]. 
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V. DISCUSSION 

Our numerical estimate of G' based on Eq.(38) sug- 
gests that a substantial enhancement of the FMR 
hnewidth, that is independent of the atomic number Z, 
should be observed in N/F/N systems. In contrast, the 
data of Mizukami et al. [8] on trilayers containing permal- 
loy films show a strong dependence of G" on Z. In fact, 
for Cm, which has the smallest Z of the A^-metals stud- 
ied there is a complete absence of a 1/d-dependent G' . 
Tserkovnyak et al. [11] propose that it is the spin ac- 
cumulation in the iV-laycr that is responsible for such 
a suppression of the ferromagnetic relaxation in copper 
layer. Note that the theory of spin-pumping assumes at 
the outset that the spin system in the iV-layer is kept 
in thermal equilibrium during the precession. For that 
one needs an efficient spin-sink mechanism. The data of 
Ref. [8] indicate that spin-lattice relaxation via spin-orbit 
coupling [21] provides the required spin sink. In fact, 
metals with larger Z exhibit generally larger measured 
values of G' . This trend is in agreement with the fact 
that the spin-lattice relaxation rate scales as Z^ [22] . 

Also the theory of G", presented in Sec. Ill, assumes 
that the electron spins in the iV-layers are in thermal 
equilibrium. This can be established either by the spin- 
lattice relaxation in the bulk, or by surface relaxation. 
One way to include these effects into the reaction field 
of Eq.(3) is to calculate the quantity (s(0,t)) using the 
Bloch equation with diffusion [23]. This equation is to 
be solved with the b.c. that allows the electron spin to 
be flipped upon collision with the surface. Such b.c. has 
been proposed by Dyson [23]. In terms of the spin density 
(s) this so called "evaporation" b.c. reads 



d{s) 2>p 
~dx ^ 4A 



is) 



(58) 



where p is the probability of spin flip to take place upon 
the reflection from the boundary, and A is the mean- 
free path in the bulk. Due to surface irregularities and 
paramagnetic surface impurities, the probability p can be 
large enough to provide the necessary spin sink even for 
layers with small bulk disorder. 

Alternatively, the spin density and the reaction field 
can be obtained from the time-dependent 2x2 matrix 
kinetic equations driven by precessing magnetization of 
the ferromagnet [24] . Such an approach is inspired by the 
work of Kambersky [25] on intrinsic damping due to spin- 
orbit coupling in bulk ferromagnets. This work invokes 
the idea of a "breathing Fermi surface": The chemical 
potential varies in response to the time-dependent per- 
turbation. However, the distribution of the electrons does 
not respond instantaneously to the perturbation. There 
is a time-lag characterized by a relaxation time r. In 
Ref. [18] we apply this idea to the case of dynamic in- 
terlayer exchange coupling. In this case it is the spacer 
electrons which are affected by the spin-dependent po- 
tential at the interfaces, and the time variation of this 
potential is due to the precession of the ferromagnetic 



moment. Thus, the relevant relaxation time is the trans- 
verse spin-relaxation time Tsp,,,, . The resulting effective 
damping field, is like Eq.(20), proportional to dM{t)/dt 
implying Gilbert damping. However, in distinction from 
Eq.(20), it is also proportional to Tgpin- 

This brings us to the question: Is there a system to 
which the ballistic theory of the present paper is applica- 
ble? We believe that the double-layer structure studied 
in Ref. [9] is a good example of such a system. Here, 
the precessing layer Fx deposits spin current into the A''- 
spacer and the second layer F2 acts as an absorber of the 
transverse component of the spin current - thus providing 
an effective spin sink. Detailed analysis of this mecha- 
nism has been presented in beautiful papers by Stiles and 
Zangwill [26,27]. These authors show that there is an os- 
cillatory, power law, decay of the transmitted transverse- 
spin current that is caused by cancellations due to a dis- 
tribution of precessional frequencies, and the rotation of 
the spin of the incoming spin upon reflection. Conse- 
quently, almost complete cancellation of the transverse 
spin takes place after propagation into the ferromagnet 
by a few lattice constants. This finding also supports our 
assumption that the excitation of transverse components 
of {s{x,t)), via s — d exchange, is confined to the N/F 
interface layer (see Appendix A). 
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APPENDIX A: DERIVATION OF EQ.(3) 

We consider a trilayer shown in Fig.l., and derive the 
reaction field from the torque equation 



Mf{t)x Hr{t) 



(Al) 



where M f{t) is the net magnetic moment of the fer- 
romagnetic film, and T{t) is the torque due to the s — d 
interaction. Since only interface regions contribute to 
this torque (see Refs.[13] and [27]), we pick a magnetic 
atom in the plane x = 0, and consider the local magnetic 
field H'^^\t) acting on its magnetic moment Ai'''-'(t). Us- 
ing Eq.(5), the expectation value of the s — d exchange 
energy of this atom is — JS^*-* (t).(s(O, t)). If we write this 
quantity as -M^'\t).H^'\t), where M^'\t) = jS^'\t), 
the local magnetic field is 



ifW(i) = ^(s(0,t)) 
7 



(A2) 
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For a square film of area _L^, the number of interface 
atoms is L? ja^. Thus, using Eq.(A2), the net torque 
contributed by both interfaces is 



2LV 
0^7 



M«(t) X (s(0,t)) 



(A3) 



Noting that the net magnetic moment of the film is 
M/ = M^^L'^d/a^, and using Eq.(A3) in (Al), yields 
Eq.(3). Similar approach has been used to deduce the 
effective field in ultrathin layers in the presence of inter- 
faces (see Eq.(1.6) in Ref.[3]). 



APPENDIX B: SPIN RELAXATION AND 
INFRARED CUTOFF 

To include spin relaxation into the theory of G', we 
start from Eq.(22) and replace the infinitesimal quan- 
tity 7] by r = h/Tg, where is the spin- relaxation time. 
Thus, the electron-hole pairs with flipped spin are as- 
sumed to relax with a frequency T/h. Similar assumption 
has been made in the theory of magnon relaxation via 
s — d interaction [28] . Moreover, we neglect the splitting 
A. Making these changes in Eqs.(22)-(24), we obtain 

jj2p roc rsi: 

ImxT{q,u^) ^ dke^^I{k,q) (Bl) 



where 

I{k,q)= [ 
Jo 



desinO 



fl^kq cos 6/m 



(B2) 



Expanding the integrand to order w, the ^-integration 
yields 



T(i. ^ 



tan X ■ 



1 -ha;2 



(B3) 



where x = fT?kq/{Tm). Using Eq.(B3) , the fc-integration 
in Eq.(Bl) yields at T = 



ImxT{q,uj) 



W?LO 

Sn'^hq 



tan '^{vpTsq) 
vpTsq 



(B4) 



1 + {vpTsqy 



Consistent with Eq.(27), this expression is equal to 

rn^uj l(\(!>i:hq) in the limit oc. The evaluation of 

the (jT-intcgral of this expression is done by approximat- 
ing tan^'^ X by -kx/2 for x <1, and by 7r/2 for x > 1. For 
r, we obtain 



dq 

—ImxT{q,(^) 



«2 



27r 



167r2fi 



lll{v FTsq2) 



(B5) 



Introducing this result into Eq.(21), the damping en- 
hancement in the presence of spin relaxation is given by 



G' 



{JMsamY q2 



87r2M 



In — 

qi 



(B6) 



where qi ^ (vpTg) ^ is the infrared cutoff mentioned at 
the end of Sec.IIIA. 
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FIG. 1: A trilayer consisting of normal metals(iV) adjacent 
to a ferromagnetic film (F) of thickness d. The s — d interac- 
tion generating the spin density s is assumed to take place in 
contact layers of thickness, o, of the order of lattice constant. 




N 



N 



X 



d 



